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This  study  is  past  of  a  continuing  effort  to  determine  the 
susceptibility  of  critical  shipboard  topside  systems  to  thermal 
radiation  and  thermal  radiation-airblast  interaction  effects  from 
nuclear  weapon  explosions.  This  report  presents  a  method  to  simply 
obtain  temperature  distribution  histories  in  system  elements  which 
are  directly  exposed  to  the  thermal  pulse  of  a  nuclear  weapon  burst. 
The  method  has  been  applied  to  four  different  elements  typically 
found  in  shipboard  systems.  Computer  programs  to  obtain  the 
temperature  response  in  these  elements  are  included. 
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INTRODUCTION 


The  susceptibility  of  shipboard  systems  to  damage  from  a 
nuclear  weapon  thermal  radiation  pulse  is  determined  by  calculating 
the  transient  temperatures  induced  in  simple  elements  of  the  system,, 
Simple  elements  chosen  are  those  which  are  the  most  vulnerable  or 
whose  failure  is  the  most  damaging  to  the  system.-  The  computation 
of  temperature  in  even  the  most  simple  elements  requires  the  use  oi‘ 
a  numerical  computer  program  based  on  finite  difference  techniques. 

Such  programs  exist  for  several  elements  such  as  circular  cylinders 
(reference  (1)),  flat  plates  (reference  (2)),  T-beams  (reference  (3)) 
and  box  beams  (reference  (4)).  However,  all  of  these  programs 
require  that  judgement  be  made  in  choosing  time  and  spacial 
increment  in  order  that  the  numerical  calculation  be  stable  and 
accurate.  Hence,  experience  in  the  use  of  numerical  computer  programs 
is  required  in  order  to  set  up  and  obtain  accurate  results  from  them. 

This  report  describes  an  efficient,  convenient  and  reasonably 
accurate  method  of  obtaining  the  temperature  response  in  selected 
elements  exposed  to  the  thermal  radiation  pulse  of  a  nuclear  weapon. 

The  method  has  been  summarized  into  simple  computer  programs  which  are 
much  less  costly  to  run  than  the  numerical  programs  and  require  no 
previous  experience  in  numerical  analysis.  The  following  steps  were 
required  to  develop  these  simple  programs.  First  a  dimensional 
analysis  is  made  of  the  element  heating  problem.  This  consists  of 
non-dimensionalizing  the  heat  conduction  equation  and  its  boundary 
conditions  which  include  heating  by  the  nuclear  weapon  thermal 
radiation  pulse  for  an  air  burst  in  the  lower  atmosphere  (reference  (6))» 
The  result  is  the  derivation  of  dimensionless  variables  which  express 
the  temperature  response  in  terms  of  a  minimum  number  of  parameters. 
Secondly,  accurate  numerical  solutions  are  carefully  obtained  from  a 
pertinent  numerical  computer  program.  Dimensionless  temperature 
distribution  histories  calculated  with  these  programs  were  punched 
on  cards  for  selected  values  of  the  dimensionless  variables  which 
were  derived  in  the  first  step.  Next,  rules  for  interpolating  or 
extrapolating  the  dimensionless  temperature  data  are  derived. 

Finally,  computer  programs  are  written  which  (1)  use  these  rules  to 
interpolate  within  the  dimensionless  temperature  arrays  (2)  convert 
the  dimensionless  temperatures  into  useful  temperature  histories  and 
(3)  print  out  all  of  the  results  in  a  convenient  format. 

This  method  of  obtaining  the  temperature  response  can,  in  principle, 
be  applied  to  any  shipboard  element.  It  will  be  illustrated  in  this 
report  by  applying  it  to  the  following  simple  elements. 
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(1)  Circular  (flinders  whose  axis  of  symmetry  is  normal  to  the 
incoming  radiation 

(2)  Plat  plates 

(3)  Circular  Cylinders  rotating  at  6  rpm  about  the  axis  of 
symmetry  normal  to  their  cross  sections 

(4)  Thermally- Thin  Elements 

The  application  to  thermally-thin  elements  requires  a  modification  of 
the  method o  Here,  the  temperature  response  is  given  by  an  energy 
balance  on  the  element,  and  the  non-dimensional  temperature  history 
is  simply  the  normalized  radiant  exposure  history  of  the  nuclear 
weapon  thermal  pulse.  Details  of  doing  this  are  given  in  a  later 
section  of  this  report. 

As  an  aid  in  deciding  if  the  solutions  generated  in  this  report 
will  be  applicable  to  a  particular  problem,  the  heat  transfer 
equations  and  boundary  conditions  for  the  four  above  elements  are 
given  in  Appendix  *  .  Also  included  in  the  appendix  is  a  figure 
(Figure  A-l)  depicting  the  geometry  of  the  elements  and  the 
positions  for  which  temperatures  can  be  found. 

A  computer  program  w«s  written  to  obtain  the  temperature 
response  of  each  of  the  above  four  elements  which  are  exposed  to 
the  thermal  radiation  pulses  of  a  nuclear  weapon.  The  use  of  these 
programs  will  be  described  a  Long  with  the  application  of  the  method 
to  these  elements .  The  computer  programs  are  not  essential  to  the 
method  of  obtaining  the  temperature  response.  For  example,  the 
dimensionless  representation  of  temperature  has  been  published  in 
graphical  form  for  flat  plates  (reference  (5))  and  for  two-dimensional 
circular  cylinders  (reference  (7)>.  It  only  remains  to  apply  the 
interpolation  rules  which  are  given  herein  to  obtain  the  desired 
temperature  histories.  But  although  these  charts  are  easy  to  use,  it 
is  time  consuming  to  compute  the  values  of  dimensionless  parameters, 
locate  the  correct  curve,  and  compute  the  temperature  histories  from 
their  dimensionless  values.  Furthermore,  it  is  frequently  necessary 
to  interpolate  between  curves  which  may  be  difficult  when  the 
interpolation  is  not  linear. 

The  techniques  developed  in  this  report  are  similar  to  previous 
work  done  at  the  University  of  California  (reference-  (8)).  The 
principal  difference  is  that  in  this  study  heating  is  caused  by  the 
nuclear  weapon  thermal  radiation  pulse  and  computer  programs  are  used 
to  retrieve  the  temperature  data. 
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SYMBOLS 


A  absorptance 

b  plate  thickness 

BCYL  basic  description  parameter  for  cylinders  (see  equation  (14)) 

BPLT  basic  description  parameter  for  plates  (see  equation  (8)) 

cp  specific  heat  at  constant  pressure 

K  thermal  conductivity 

H  irradiance 

Hm  maximum  irradiance 

Q  radiant  exposure 

Qr  absorbed  radiant  exposure 

QT  total  radiant  exposure 

r  cylinder  radius 

r*  dimensionless  cylinder  radius  (r/Ro) 

Ri  cylinder  inner  radius 

Ri*  dimensionless  cylinder  inner  radius  (Ri/Ro) 

Ro  cylinder  outer  radius 

t  time 

tm  time  at  which  maximum  irradiance  of  a  nuclaar  weapon  thermal 
pulse  occurs 

t*  dimensionless  time  (t/tm) 

T  temperature 

T  average  temperature 

T*  dimensionless  temperature  (see  equations  (5)  and  (9)) 

TI  initial  temperature 

X  distance  into  plate,  measured  from  its  surface 

X*  dimensionless  distance  into  plate  (x/b) 

a  thermal  diffusivity  (K/^cp) 

u  angular  velocity 

p  density 

9  included  cylinder  angle,  measured  from  the  most  forward  point 
on  the  cylinder 

9*  dimensionless  cylinder  angle  ($/*) 
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THE  DIMENSIONLESS  REPRESENTATION  OF  THERMAL  RESPONSE 


A  dimensionless  analysis  of  the  heat  conduction  equation  with 
thermal  radiation  heat  input  from  a  nuclear  weapon  detonation  was 
made  for  a  plate  (reference  (5))  and  for  a  circular  cylinder 
(reference  (7))„  The  resulting  dimensionless  groups  will  not  be 
derived  again  here,  instead  they  will  be  simply  stated  along  with 
the  assumptions  required  to  produce  them.  The  thermal  radiation 
heating  is  given  by  the  normalized  nuclear  weapon  thermal  radiation 
pulse  given  in  reference  (6)  and  shown  in  Figure  1.  This  curve 
represents  the  heat  energy  directed  to  any  body  which  is  directly 
exposed  to  a  nuclear  weapon  burst.  The  total  heating  per  unit 
area  (the  radiant  exposure)  received  by  the  body  up  to  time,  t,  is 
simply  the  integral  of  the  thermal  pulse  up  to  that  time,  i.e.. 


Q  * 


H  (t)dt. 


(1) 


Since  the  weapon  thermal  pulse  is  given  in  normalized  form,  the 
radiant  exposure  is  more  properly  expressed  as: 


Where  QT  =  2.6  Hmtm  is  the  total  radiant  exposure 


•f?%  1  ■*  &**  ynv«^f»»^<  gy^?C^gp5ggy3S"g^r  iggggy VB**EFlQf&B%!S&QfV!!t^£r&Gem&& ffggg^gjg^g|^}e ^ 
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Only  that  portion  of  the  radiant  exposure  which  is  actually 
absorbed  by  the  target  element  will  cause  its  temperature  to  rise. 

The  quantity  of  heat  that  is  absorbed  depends  on  the  absorptance  of 
the  target  surface  material  and  on  the  radiative  "shape  factor" 
which  describes  how  the  target  is  oriented  with  respect  to  the  thermal 
radiation  source.  I*-  was  assumed  that  the  absorptance  is  a  constant 
for  any  given  target  material  and  that  the  target  elements  are 
directly  exposed  to  the  weapon  explosion.  Hence,  the  amount  of 
thermal  energy  absorbed  is 

Qr  =  A  Q  for  a  plate  (3) 

and 

Qr  =  A  Q  cos  e  for  a  circular  cylinder  (4) 

Equation  (3)  or  (4)  represents  the  boundary  condition  for  radiant 
heat  received  at  the  outer  surfaces.  Other  boundary  assumptions  are 
that  no  heat  is  received  at  the  inner  surfaces  and  no  convective 
or  radiative  cooling  occurs  at  either  surface.  It  was  also  assumed 
that  the  thermophysical  properties  (A,  K,  p  ,  cp  and  a )  do  not  depend 
on  the  temperature  of  the  element.  It  was  shown  in  reference  (7) 
that  these  assumptions  have  little  effect  on  the  temperature  history 
of  shipboard  elements  in  practical  applications.  Using  these  assump¬ 
tions  and  boundary  conditions,  dimensionless  variables  were  chosen 
such  that  a  minimum  number  of  parameters  were  required  to  define  the 
temperature  distribution  history.  The  dimensionless  groups  found 
were  the  following. 


(a)  For  plates 


(1)  Dimensionless  Temperature  T*  =  (T-TI) / (AHmb/K)  (5) 

(2)  Dimensionless  Time  t*  =  t/tm  (6) 

(3)  Dimensionless  Distance  into  Plate 

X*  =  X/b  (7) 

(4)  Basic  Description  Parameter  „ 

BPLT  =  acm/'c  (8) 


The  temperature  rise  of  plates  was  given  in  reference  (5 )  by 
plotting  graphs  of  dimensionless  temperature  versus  dimensionless 
time  for  12  values  of  the  basic  description  parameter.  Figures 
were  presented  for  three  values  of  dimensionless  distance  which 
correspond  to  the  surface,  midplane  and  backface  of  the  plate. 

(fc)  For  circular  cylinders 


(1)  Dimensionless  Temperature  T*  =  (T-TI)/ (AHmRo/K)  (9) 

(2)  Dimensionless  Time  t*  =  t/tm  (10) 

(3)  Dimensionless  Radius  into  Cylinder 

r*  =  r/Ro  (11) 

(4)  Dimensionless  Inner  Radius  * 

Pi  =  Ri/Ro  (12) 

(5)  Dimensionless  Cylinder  Angle 

»*  -  O/ir  (13) 

(6)  Basic  Description  Parameter 

BCYL  =  atm/Rc2  (14) 
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The  temperature  rise  of  circular  cylinders  was  given  in  reference  (7) 
by  plotting  graphs  of  dimensionless  temperature  versus  dimensionless 
time  for  ten  values  of  the  basic  description  parameter.  Figures 
were  presented  for  three  values  of  the  dimensionless  radius,  five 
values  of  the  dimensionless  inner  radius,  and  six  values  of  the 
dimensionless  cylinder  angle  (a  total  of  90  figures). 

A  dimensional  analysis  similar  to  that  made  for  circular 
cylinders  could  be  made  for  rotating  circular  cylinders  and  an 
additional  parameter (s)  derived  to  account  for  rotation.  However, 
this  will  not  be  done  since  the  temperature  histories  of  this  report 
were  obtained  for  one  rotational  speed  only  (6  rpm) .  Therefore,  the 
dimensionless  representation  of  the  stationary  cylinders  will  be  used 
to  describe  the  single-speed  rotating  cylinders  discussed  herein. 


INTERPOLATION  OF  THE  DIMENSIONLESS  DATA 

Dimensionless  temperature  histories  have  been  computed  for 
plates,  circular  cylinders  and  circular  cylinders  rotating  at  6  rpm 
which  are  directly  exposed  to  the  thermal  radiation  pulse  of  a 
nuclear  weapon.  The  calculations  were  made  using  the  previously 
mentioned  computer  programs  which  are  based  on  finite  difference 
techniques.  Accuracy  of  the  calculations  was  insured  by  a  careful 
choice  of  time  and  distance  increments  to  be  used  in  the 
computations.  (See  reference  (7)  for  a  discussion  of  the  accuracy 
of  numerical  solutions.)  A  large  number  of  computer  runs  were  made 
and  the  temperature  histories  computed  were  converted  to  dimensionless 
values  and  punched  on  computer  cards.  Cards  were  punched  for 
selected  values  of  the  dimensionless  variables.  These  selected 
values  are  given  in  Table  1  for  plates,  circular  cylinders  and 
rotating  circular  cylinders.  The  dimensionless  temperature  data  was 
organized  such  that  the  dimensionless  temperatures  are  given  in  terms 
of  the  dimensionless  time  with  all  of  the  remaining  dimensionless 
par  *ime ters  varied  in  turn. 

Since  the  parameter  values  of  a  given  problem  may  not  coincide 
with  those  for  which  a  calculated  solution  exists,  it  is  generally 
necessary  to  interpolate  or  extrapolate  among  those  parameter  values 
for  which  calculations  have  been  made.  The  procedure  was  to  find 
the  dimensionless  temperature  history  for  only  the  values  of 
dimensionless  time  (t*)  and  dimensionless  geometric  variables  (r*,  «* 
or  X*)  which  are  listed  in  Table  1.  Hence  it  is  necessary  to 
interpolate  among  the  remaining  dimensionless  variables  only 
(BCYL  and  Ri*for  the  cylinders  and  BPLT  for  the  plate). 

Rules  for  the  interpolation  of  dimensionless  temperature  data  for 
cylinders  were  derived  in  reference  (7).  A  distinction  is  made 
when  the  element  is  thermally-thick  or  not.  An  element  is  defined 
to  be  thermally-thick  at  any  time  if  its  inner  surface  temperature 
rise  is  less  than  one  percent  of  its  outer  (exposed)  surface 
temperature  rise  at  that  time.  The  following  rules  will  be  used  for  the 
interpolation  or  extrapolation  of  dimensionless  temperature  data. 
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(1)  Thermally-thick  elements 

(a)  Cylinders  T*~<jBCYL  .  t*  (15) 

(b)  Plates  .  t*  •  (16) 

-  ■  '  6 

(2)  Non  thermally-thick  elements  ’ 

(a)  Cylinders  T*^T^7rI*  (17) 

(b)  Plates  ;  T*~  BPLT  -  t*  .  (18) 

i 

The  above  equations  show  that  linear  interpolation  is  used  only 
for  the  case  of  non  thermal lyr thick  plates  (see  equation  (18) ) .  In 
dealing  with  thermally- thick  elements#  it  is  necessary  to  use  square 
root  interpolation  on  the  values  of  the  basic  description  parameter. 

The  non  thermally- thick  cylinders  are  the  most  difficult  to  deal  with ; 
since  double  interpolation  on  the  parameters  BCYL  and  Ri*is  required 
as  specified  by  equation  (17).  Since 'equations  (15)  through  '  f  18)  are 
only  approximate  representations,  of  the  behavior  of  the  dimensionless 
temperature  histories,  the  interpolation  rules  will  ■ yield, only  an 
approximate  result.  The  expected  accuracy  in  interpolating  the 
dimensionless  data  is  given  in  a  lateir  section  of  this  report  to  be 
less  than  +5  percent  for  the  data  banks  used  in  the  computer  programs 
of  this  report.  i 

DESCRIPTION  AND  USE  OP  THE  COMPUTER  PROGRAMS 

i  * 

I  5 

The  computer  programs  used  ,to  retrieve  the  temperature 
distribution  histories  require  as  input  the  material  properties  and 
dimensions  of  the  element  along  with  the  total  radiant  exposure 
and  time  at  which  the  maximum  irradiance  occurs  for  the  weapon 
yield  of  interest.  Also  input  is  the  dimensionless  temperature  data 
which  is  in  the  form  of  decks  of  punched  cards  or  a  magnetic  tape. 

From  this  input,  the  dimensionless  parameters  are  computed  and 
interpolation  or  extrapolation  of  the  dimensionless  temperature  data 
is  accomplished  using  the  previously  discussed  interpolation  rules 
(see  equations  (15)  through  (18) ) .  The  dimensionless 1  temperatures 
calculated  are  then  made  dimensional  and  subsequently  printed  out 
in  a  convenient  format.  Four  computer  programs  have  been  written, 
one  for  stationary  circular  cylinders,  one  for  plates,  one  for  circular 
cylinders  rotating  at  6  rpm  and  one  for  thermally- thin  elements.  These 
are  listed  in  Figures  2  through  5,  Respectively.  Fortran  decks  of 
these  programs  and  punched  cards  or  magnetic , tape  of  the  dibensionless 
temperature  data  may  be  obtained  from  NOL.  The  definition  of  the 
variables  that  appear  in  the  input  of  these  programs  and  their  units 
ire  given  in  Appendix  B.  Note  that  all  of  these  variables  have  the 
same  definition  regardless  of  in  which  program  they  are  used.  The  1 
following  is  a  description  of  each  program  in  turn,  including 
instructions  on  punching  the  input  data  cards  required  to  run  the 
program. 

(1)  Retrieval  of  temperature  histories  in  stationary  circular  . 
cylinders  (Program  CYLDAT) 
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Only  two  or  three  data  input  cards  are  required  to  run  this 
program.  All  of  the  stationary  circular  cylinder  dimensionless 
temperature  data  has  been  placed  on  a  magnetic  tape  which  makes  their 
retrieval  orach  simpler.  The  variables  and  their  format  required  to 
-‘run  the  program  are  given  in  the  following  table 

CARD  VARIABLE  FORMAT 

1  *  XK,  RHO,  CP,  TM,  RI,  RO,  TI,  H  7E10.4,  110 

2.3  QT,  ABSR,  TIME (m) ,  M=  1,H  2E10.4,  12F5.2/12F5.2 

3.4  ,  Blank  if  no  further  runs,  or  repeat  cards  1-3  for  a  new  run 

Figure  2  is  a  complete  Fortran  listing  of  this  program.  The  choice 
of  either  two  or  three  data  input  cards  arises  from  the  choice  of  the 
number  (N)  of  dimensionless  time  values  (TIME)  desired.  Note  that 
up>  to -24  values  of  dimensionless  time  may  be  chosen  from  the 
allowable  values  given  in  Table  1.  If  12  or  fewer  are  chosen,  they 
will  be  punched  on  card  2  and  card  3  is  not  present.  A  feature  of 
this  program  is  that  if  all  24  values  of  dimensionless  time  are 
'  ’sired,  the  variable  N  is  made  zero  and  no  values  of  dimensionless 
time. need  be  specified  on  card  2.  (Card  3  is  omitted). 

(2)  Retrieval  of  temperature  histories  in  flat  plates  (Program 
PLTlD) . 

The  initial  data  input  cards  required  to  run  this  program 
are  very  similar  to  those  of  the  previous  program.  A  major 
difference  in  input  arises  because  the  dimensionless  temperature  data 
for  plates  is  available  on  punched  cards  only.  Hence,  card  decks  of 
dimensionless  temperature  data  roust  be  read  in  before  the  initial 
two  pr  three  data  input  cards.  The  dimensionless  temperature  data 
consists'  of  12  decks  of  punched  cards  (15  cards  per  deck)  which 
correspond  to  the  12  selected  values  of  the  parameter  BPLT  (see 
:  Table  1).  In  each  deck,  the  dimensionless  temperatures  are  punched 
for  the  25*  selected  values  of  dimensionless  time  and  the  five 
selected  values  of  dimensionless  distance  into  the  plate  (see 
Table  1).  The  input  data  required  to  run  this  program,  including 
the  variables  and  their  format  for  the  initial  data  cards,  are  given 

in  the  following  table. 

» 

CARD  VARIABLES  FORMAT 

1  XK,  RHO,  CP,  TM,  THICK,  TI,  N  6E10.4,  110 

2,3;  QT,  ABSR,  TIME  (M),  M  =  1,  N  2E10.4,  12F5.2/13F5.2 

Blank  card  if  no  further  runs,  or  repeat  above  cards  for  a 
new  run 


*  One  ektra  dimensionless  temperature  corresponding  to  t*  =  0  was 
included  in  the  plate  data. 
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Figure  3  is  a  complete  Fortran  listing  of  this  program.  The 
variable,  H,  nay  also  be  ndt  zero  to  obtain  temperatures  for  all 
of  the  allowable  values  of  tine  as  described  for  the  previous  program. 

(3)  Retrieval  of  temperature  histories  in  circular  cylinders 
rotating  at  6  rpn  (program  ROTCYL) 

the  use  of  this  progran  is  similar  to  that  of  the 
preceding  program  in  that  it  requires  two  or  three  initial 
data  cards  followed  by  a  deck(s)  of  the  dimensionless  temperature 
data.  Here,  there  are  also  12  possible  decks  of  dimensionless 
tenperature  data  corresponding  to  the  three  selected  values  of  the 
parameter,  at,  and  four  selected  values  of  the  parameter  BCYL  (see 
Table  1).  In  each  deck,  the  dimensionless  temperatures  are  punched 
for  the  24  selected  values  of  dimensionless  time,  the  five  selected 
values  of  dimensionless  radius  and  the  six  selected  values  of  dimensionle 
cylinder  angle  (see  Table  1) .  To  set-up  this  program,  it  is  first 
necessary  to  calculate  the  values  of  Rx*and  BCYL  for  the  problem 
under  consideration.  These  values  must  be  inspected  to  determine 
if  interpolation  or  extrapolation  of  the  allowable  values  of  Rx*or 
BCYL  is  required.  A  variable,  KK,  is  input  to  indicate  vhich 
variables  are  to  be  interpolated  or  extrapolated.  KK  is  chosen 
from  the  following  table. 

KK  JfEJUlIHG 

1  Ho  interpolation  or  extrapolation  of  either  parameter 

2  Ho  interpolation  or  extrapolation  cm  Rif  interpolate  on  BCYL 

3  HO  interpolation  or  extrapolation  on  BCYL,  interpolate  or 

extrapolate  on  Ri* 

4  Interpolate  on  BCYL,  interpolate  on  Ri* 

5  Interpolate  on  Ri*,  extrapolate  cm  BCYL 

6  Ho  interpolation  or  extrapolation  on  BCYL,  extrapolate  cm  Ri* 

7  Interpolate  on  BCYL,  extrapolate  on  Ri* 

The  program  for  stationary  cylinders,  previously  described,  also  does 
all  of  these  extrapolations  or  interpolations  (except  number  5)  but 
there  the  correct  value  of  KK  is  selected  by  the  program  itself.  Two 
variables,  I  and  J,  are  input  to  indicate  which  dimensionless 
tenperature  decks  are  to  be  used  by  the  program.  I  indicates  the 
values  of  BCYL  and  is  chosen  from  the  following  table. 

BCYL  .010  .025  .050  .100 

I  12  3  4 

When  interpolation  is  required,  I  will  correspond  to  the  larger 
value  of  BCYL.  If  extrapolation  beyond  BCYL  =  .100  is  required,  I 
is  set  equal  to  four,  but  extrapolation  below  BCYL  =  .01  is  not 
allowed.** 


**  Additional  dimensionless  data  decks  could  be  generated  for  this 
case  from  the  numerical  programs  previously  mentioned. 
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Hie  integer  J  indicates  the  value  of  Ri*  and  is  chosen  froe  the 
following  table. 


Ri*  0.95  0.90  0.80 

J  1  2  3 

when  interpolation  is  required,  J  will  correspond  to  the  analler 
value  of  Ri*.  If  extrapolation  below  Ri*  =  0.8  is  required,  J  is 
set  equal  to  three,  but  extrapolation  above  Ri*  =  0.95  is  not 
allowed.  The  total  input  data  required  to  ran  this  program, 
including  the  variables  and  their  foreat  for  the  initial  data  cards, 
are  given  in  the  following  table. 

CARO  VARIABLE  FORMAT 

1  XK,  RH0,  CP,  TM,  RI,  R0,  TI,  KK,  I,  J,  K  7E10.4,  2(13,12) 

2,3  QT,  ABSR,  TIME  (M),  1*0.,  H  2E10.4,  12F5.2/12F5.2 

1-4  decks  of  dimensionless  temperature  data  (120  cards  per 
deck).  Blank  card  if  no  further  runs,  or  repeat  the  above  cards 
for  a  new  ran. 

The  d  imens  ion  less  temperature  decks  are  loaded  beginning  with  the 
deck  having  the  smallest  value  of  I  (smallest  BCYL  value)  and 
largest  value  of  J  (smallest  Ri*  value).  Succeeding  decks  will  have 
the  smaller  value  of  J  and  larger  value  of  I  followed  by  the  smaller 
I  and  smaller  J  and  finally  by  the  larger  I  and  smaller  J  (omit 
any  of  these  not  required) .  Figure  4  is  a  complete  Fortran 
listing  of  this  program.  The  variable,  H,  may  also  be  made  zero  to 
obtain  temperatures  for  all  of  the  selected  values  of  time  as  was 
done  in  the  two  previous  programs. 

(4)  Retrieval  of  temperature  histories  in  therraally-thin 
elements  (Program  “THNEL") 

This  program  finds  the  average  teaperature  history  of  a 
thermally- thin  element,  or  more  precisely,  of  an  element  in  which 
no  heat  is  internally  conducted.  Target  elements  are  usually  defined 
to  be  thermally- thin  if  there  is  less  than  a  ten  percent  difference 
between  the  front  and  rear  surface  temperatures  and  thermally- thick 
if  the  rear  surface  temperature  is  less  than  one  percent  of  the 
surface  temperature .  Thus,  in  order  to  use  this  program,  it  is 
necessary  to  know  where  internal  heat  conduction  effects  are 
important.  Figure  6  shows  the  thermally- thin  and  the  thermally- thick 
regions  for  a  plate  in  terms  of  the  basic  description  parameter,  BPLT, 
and  the  dimensionless  time,  t*.  For  plates,  the  temperature 
histories  retrieved  by  this  program  sure  average  temperature 
histories  but  these  are  an  accurate  representation  of  the  temperature 
history  in  the  thin  plate  region  only. 

Figure  7  shows  the  thermally- thin  and  thermally-thick  regions 
of  a  stationary  circular  cylinder  in  terms  of  the  basic  description 
parameter,  BCYL,  and  both  the  dimensionless  time,  t*,  and  the 
dimensionless  inner  radius,  Ri*.  Also  indicated  on  Figure  7  are 
regions  where  angular  (circumferential)  conduction  is  important. 
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For  cylinders,  the  twprratntt  histories  retrieved  by  this  program 
are  the  average  temperature  only  in  the  region  where  angular 
conduct  ion  is  negligible.  Thus,  these  temperature  histories  are 
accurate  representations  of  the  temperature  history  in  a  cylinder 
only  in  thse  regions  abew  in  Figme  7  den  there  is  no  angular 
conduction  and  the  cylinder  is  tfokmally-thia.  In  this  case,  the 
teaperature  history  at  any  point  on  the  cylinder  is  the  retrieved 
value  Multiplied  by  cos  f  where  •  is  the  cylinder  angle  neasured 
tram  the  nost  forward  point  on  the  cylinder. 

The  average  teaperature  history  is  found  by  neans  of  a 
siaple  energy  balance,  i.e.. 


T  *■ 


AQ 

»  cp  b 


+  n 


(19) 


The  radiant  exposure,  Q,  was  found  fay  integrating  equation  (2)  and 
normalizing  the  result  with  the  total  radiant  exposure,  QT.  Since 
Q/QF  in  tens  of  t*  is  a  constant  fraction  four  nuclear  weapons  of 
all  yields  (i.e.,  all  values  of  ta),  it  was  stored  as  peraanent 
data  in  the  program.  This  fraction  is  subsequently  used  in  the 
program  to  compute  the  radiant  exposure  and  the  teaperature  histories 
an  calculated  by  use  of  equation  (19). 

(My  one  data  input  card  is  required  to  run  this  program.  The 
variables  and  their  format  for  this  card  are  given  in  the  following 
table. 

CARD  VARIABLES  FORMAT 

1  RHO,  CP,  TM.  THICK,  TI,  QT,  ABSR  7B10.4 

2  Blank  card  if  no  further  runs,  or  repeat  card  1  for  a  new  run 

Figure  4  is  a  complete  Fortran  listing  of  the  program.  This  program 
provides  teaperature  histories  for  all  26  selected  values  of 
dimensionless  time  given  in  Table  1.  The  program  has  been  simply 
modified  to  compute  the  teaperature  histories  in  rotating  cylinders. 
The  modification  consists  of  replacing  the  Q/QT  values  with  values 
obtained  by  integrating  equation  (2)  with  H  replaced  by  H  cos  »  . 

Here  t  is  the  time- varying  angle  between  the  incoming  radiation  and 
a  line  normal  to  the  cylinder  at  its  most  forward  point  as  it 
rotates.  Figure  7  for  stationary  cylinders  can  be  used  as  a  guide 
in  finding  the  region  where  a  rotating  cylinder  is  both  thermally-thin 
and  has  negligible  angular  conduction.  In  Figure  7,  the  thin  cylinder 
region  increases  but  the  negligible  angular  conduction  region 
decreases  as  the  rotation  rate  is  increased. 


Appendix  C  presents  examples  for  the  set-up  and  use  of  each 
of  the  computer  programs  described  in  this  section. 
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ACCURACY  OP  THE  MBTBtJD 


The  present  Method  of  obtaining  temperature  distribution 
histories  in  target  elements  consists  primarily  of  solving  some 
representative  problems,  storing  these  solutions  in  dimensionless 
forms  and  retrieving  the  solution  for  a  given  problem  by  interpolating 
the  dimensionless  data.  The  calculations  for  the  representative 
problems  were  made  by  the  large  numerical  computer  programs  based  on 
finite-difference  techniques.  These  are  accurate  solutions  because 
spacial  and  temporal  increments  mere  carefully  selected.  These 
calculations  considered  elements  of  constant  theemopb  ■  sical  properties 
(K,  0  and  cp)  and  neglected  any  radiation  or  convective  cooling  cm  the 
inner  or  outer  surfaces.  However,  these  restrictions  do  not  seriously 
limit  the  utility  of  the  present  method,  primarily  because  they  have 
a  negligible  effect  for  temperature  rises  which  occur  in  shipboard 
applications  (see  reference  (7)).  Por  example,  aluminum  alloys 
vhich  are  used  in  many  structures  cannot  be  heated  to  more  than  a  few 
hundred  degrees  centigrade  or  they  will  rapidly  lose  much  of  their 
strength.  Also,  the  error  involved  in  using  constant  thermophysical 
properties  can  be  reduced  by  selecting  property  values  midway  between 
the  initial  and  expected  maximum  temperature  values. 

Hence,  there  remains  only  the  error  introduced  by  the 
interpolation  rules  (equations  (15)  -  (18)).  This  error  fan  only  be 
estimated  because  the  interpolation  rules  are  only  approximate 
representations  of  dimensionless  temperature  behavior.  Figure  8 
shows  this  error  by  comparing  the  results  of  an  accurate  numerical 
calculation  of  dimensionless  temperature  history  with  the  result 
obtained  by  '-he  interpolation  rules.  Also  shown  is  the  result  of  a 
linear  interpolation  of  the  dimensionless  data.  In  this  example,  a 
circular  cylinder  having  dimensionless  parameter  BCYL  and  Ri*  values 
of  0.075  and  0.85,  mas  analyzed.  These  values  were  chosen  to 
represent  a  "worst  case"  for  interpolation  because;  (1)  interpolation 
on  both  dimensionless  parameters  is  required  and  (2)  there  are 
relatively  large  differences  in  dimensionless  temperature  values  between 
the  selected  values  of  the  dimensionless  parameters  which  are  used 
for  interpolation.  Figure  8  shows  that  the  interpolated  result  is 
within  approximately  two  percent  of  the  numerical  calculation.  Also 
the  straight- forward  linear  interpolation  is  within  approximately 
five  percent  of  the  numerical  result.  The  reason  that  even  the  linear 
interpolation  provides  such  good  results  is  that  the  dimensionless 
data  provided  is  sufficiently  plentiful  that  only  a  small  error 
is  possible  in  the  interpolations.  Based  on  this  example  and  other 
examples  given  in  reference  (7),  it  is  concluded  that  the  average 
deviation  in  temperature  history  in  using  the  present  interpolation 
method  is  less  than  +  5  percent  when  compared  with  an  accurate 
numerical  calculation . 
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SUMMARY  AMD  CCHCMJSICHS 

this  report  has  presented  a  method  of  easily  obtaining  the 
tempera  tore  distribution  history  in  some  simple  elements  iAich 
have  been  directly  exposed  to  the  thermal  radiation  pulse  of  a 
nuclear  weapon.  Although  this  method  requires  the  use  of  computer 
programs,  these  programs  are  much  easier  to  set-up  and  run  than 
the  numerical  programs  which  compete  temperature  distribution 
histories.  Also,  these  simpler  programs  do  not  require  that  judg¬ 
ments  be  made  in  choosing  input  variable  values.  Furthermore,  the 
simpler  programs  always  require  much  less  computer  time  to  produce 
answers.  (Typical  run  time  is  less  than  ten  seconds  per  problem  on 
the  GDC  6400  computer. } 

The  main  restrictions  on  the  present  method  are  that  convective 
or  radiative  cooling  is  not  accounted  for  and  that  the  thermopfaysical 
properties  are  constant.  These  restrictions  do  not  seriously  effect 
the  results  of  this  method  provided  that  the  temperature  rises  generated 
are  limited  to  a  few  hundred  degrees  centigrade.  Since  nary 
shipboard  structural  elements  are  constructed  of  aluminum  alloys 
which  are  only  useful  in  this  temperature  range,  the  restrictions  do 
not  seriously  limit  the  application  of  the  method.  The  present  cooputex 
programs  are  restricted  to  the  range  of  dimensicmle  is  parameters  shown 
in  Table  1.  This  is  not  a  serious  limitation  because  a  wide  range 
of  weapon  yields,  radiant  exposures  and  element  sizes  are  included 
in  the  given  range  of  dimensionless  parameters.  This  limitation 
can  be  remedied  by  providing  additional  dimensionless  temperature 
data  for  a  wider  range  of  dimensionless  parameters .  Also,  extrapolation 
of  the  data  is  allowed  but  must  be  performed  with  care  to  avoid  error. 
The  method  of  this  report  therefore  provides  a  reliable,  inexpensive 
and  reasonably  accurate  way  of  obtaining  the  temperature  distribution 
history  of  many  practical  shipboard  structural  elements  which  have 
been  directly  exposed  to  the  thermal  radiation  pulse  of  a  nuclear 
weapon  detonation. 
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PROGRAM  CYLBAT ( INPUT ,TAPEl , OUTPUT ) 

DIMENSION  X(10),Y(5),W(8),Z(24),R(5),C(24,8,5),D(24,8,5),E(24,8,5) 
1  ,TIME  (24 1 , TEMPI  (5)  ,TEMP2  (5)  ,F  (24,8,5) 

DATA  X/.G001, .00025, .0005, .001, .002, .005, .010, .025, .050, .100/ 

DATA  Y/0. 9, 0.8, 0.6, 0.4, 0.0/ 

DATA  W/0. 0,15. 0,30. 0,45. 0,60. 0,90. 0,135. 0,180,0/ 

DATA  Z/0. 5,0. 75, 1.0, 1.2, 1.4, 1.6, 1.8, 2. 0,2. 2, 2. 5, 2. 8, 3. 25, 4. 00, 5.0, 
+6.0,7.0.8.0,9.0,10.0,11.0,12.0,14.0,16.0,20.0/ 

1  READ  1001,XK,RHO,CP,TM,RI,RO,TI,N 

1001  FORMAT (7P10. 4, 110) 

IF (NoEQ.O)GO  TO  2 

11  IF  (XK)  120,120,9 
9  READ  1002, QT,ABSR,  (TIME  (M)  ,M=1,N) 

1002  FORMAT  (2E12 .4, 10F5 . 2/15F5 . 2 ) 

502  FORMAT (1H0,2X,56H  DIMENSIONLESS  TEMPERATURES  FOR  A  DIMENSIONLESS  T 
1IME  OF  ,F5.2,30X,6H  TIME= , F6  <  3 , 7HSECONDS ) 

503  FORMAT (1H0,11H,  CYL  ANGLE*,  F4 . 0 , 3HDEG, 3X, 5HT*  =  , 5F7 .4, 5X, 10HCYL  AN 
1GLE=,F4.0.3HDEG,2X,5HT*  =,5F7.4) 

504  FORMAT (1H  ,21X,5HIV/RO=,5F7.4,24X,5HR/RO=,5F7.4) 

505  FORMAT (1H  ,21X, 5HT(F)=, 5F7.1, 24X, 5HT (F)=, 5F7 .1) 

GO  TO  4 

2  IF (XK)  120,120,7 

7  READ  1002 ,  QT , ABSR 
DO  6  M=l, 24 
6  TIME(M)=Z(M) 

Nl=24 

4  ALPHA=XK/ (RHO*CP) 

A=ALPHA*TM/ )  RO*RO ) 

B=RI/R0 
PRINT  401, A, B 

401  FORMAT (1H1,31H  BASIC  DIMENSIONLESS  PARAMETERS, 20X,  5HBCYL=,F7 .6, 
15X,6HRI/RO=,F4.3) 

PRINT  402,XK,RHO,CP 

402  FORMAT (1H,  28H  MATERIAL  THERMAL  PROPERTIES, 4Xf 2 1HTHERMAL  CONDUCTIV 
1ITY= , F6 . 3 , 12HCAL/ SEC-CM-C , 4X, 8HDENSITY* , F5 . 2 , 7HGMS/CM3 , 4X, 14HSPECI 
2FIC  HEAT= , F4 . 3 , 8HCAL/GM-C ) 

PRINT  403,TM,QT 

403  FORMAT (1H  , 31H  WEAPON  THERMAL  CHARACTERISTICS, 10X, 2 3HTIME  TO  MAX  I 
1RADDIANCE* , F7 . 4 , 7HSEC0NDS , 9X, 2 3HTOTAL  RADIANT  EXPOSURE* , F8 „ 2 , 7HCAL 
2/ CM2 ) 

PRINT  404,RI,RO,TI 

404  FORMAT ( 1H  , 27H, 13HINNER  RADIUS* , F6 . 3, 11HCENTIMETERS , 4X, 13HOUTER  RA 
1DIUS=,F6. 3, 11HCENTIMETERS , 4X, 20HINITIAL  TEMPERATURE*, F5.1, 5HDEG-C) 

HM=QT/(2.60*TM) 

DEM=ABSR*HM*RO/XIC 
R  (1)=B 

S=(1.0-B)/4.0 
DO  3  K=2 , 5 
Kl=K-l 

3  R (K)*R)Kl)+S 
KK=0 

Figure  2.  Computer  Program  to  Obtain  the  Temperature  History  Distri¬ 
bution  in  Circular  Cylinders 
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TOLXI*0. 00009 
T0LXE-0.il 
TOLY=0  o 905 

IP (A.LE.TOIXI .OR.A.GE.TOLXE)KK*KK+l 
IP  (B.GE.TOLY)KK*KK+2 
IF(KK.GE.l)KK*KK+3 
5  DO  10  1*1,10 
T0L1*0.9*X(I) 

IP  (A-TOLl)  15,15,8 
8  TOL2«lol*X(I) 

IF (A-TOL2 )  15,15,10 
10  CONTINUE 
15  DO  20  0*1,5 
TOLl*0.995*Y (J) 

IP (B-T0L1)  18,24,24 
18  TOL2*1.005*Y (J) 

IF (B-TOL2)  20,24,24 
20  CONTINUE 

24  IF(KK)  25,25,40 

25  TOLXI=0.9*X(I) 

TOLXE»l.l*X(I) 

IF ( A . LE . TOLXI , OR . A . GE . TOLXE ) KK*KK+ L 
TOLY 1*0 „ 996  *Y ( J ) -0  o  001 
TOLYE-1. 004 *Y(J) +0.001 
IF ( B . LE . TOLY I .OR.B.GE. TOLYE ) KK*  *2;  2 
40  KK*KK+1 

IF (KK.LEo4)GO  TO  45 
J»2 

KK*KK+1 

DO  46  KJ*1, 10 

TOL1=0.9*X(KJ) 

TOL2=l.l*X(KJ) 

IF ( A  o  GE . TOLI . AND . A . LE . T OL2 ) KK=KK- 1 
46  CONTINUE 

45  GO  TO  950, 55, 70, 60, 70, 70, 60 )KK 
50  N*10*(J-1)+I 
DO  52  L=1,N 

52  READ (1) ( ( (E (Kl,K2,K3) , Kl=l, 24) ,K2*1, 8) , K3=l, 5) 
GO  TO  88 

55  N»10*(J-1)+I 
NN=N~1 

DO  56  L=1,NN 

56  READ (1) ( ( (C (Kl,K2,K3) ,Kl*l, 24) ,K2*1, 8) ,K3=1, 5) 
READ(l) (( ((D(K1,K2,K3),K1=1,24),K2=1,8),K3=1,5) 
GO  TO  88 

60  N=10*  (J-2 )  +  l 
NN=N-1 
DO  65  L*1,NN 

65  READ ( 1 )  ( ( (C (Kl,K?,X3) ,Kl=l, 24) ,K2=1, 8) ,K3=1,  5) 
READ(l) ( ( (D(K1,K2,K3) ,Kl»l,24) ,K2=1,8) ,K3=1,5) 
GO  TO  80 
70  N=10*  (J-2)  +  l 
DO  75  L=1,N 
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75  READ(l)  (<(C(Kl,K2,K3),Kl=l,24),K2=l,8),K3=l,5> 

HK=S+1 
K=N+10 
DO  77  L=NH,N 

77  READ(l) (((E'KlrK2,K3),Kl=l,24),K2=l,8),K3=l,5) 

GO  TO  88 
80  BJK-N+1 
H=S+9 

DO  85  L=H1I,H 

85  READ(l) (((E(Kl,K2,K3),Kl=l,24),K2=l,8),K3=l,5) 

READ(l) ( ( (F (Kl,K2,K3) ,Kl=l,24) ,K2=1,8) ,K3=1,5) 

88  11=1 

DO  110  L=l,24 

IF (TIME (II) -Z (L) )  90,90,110 

90  GO  TO  (108,  91,98, 98,  91, 98, 98) KK 

91  IF(C(L,1,1).GE.0.01*C(L,1,5))  GO  TO  92 

FRAC= (SQRT (A) -SQRT (X (1-1) ) )/ (SQRT (X (I) ) -X (1-1) ) ) 

GC  TO  93 

92  FRAC=  (A— X(I— 1) )/  (X  (I)-X  (1-1) ) 

93  J1=L 

DO  94  J2=l,8 
DO  94  J3=l,  5 

94  E(J1,J2,J3)=C(J1,J2,J3)+FRAC*(D(J1,J2,J3)-C(J1,J2,J3) ) 
GO  TO  108 

98  IF(C(L,1,1),LE#0.01.*C(L,1,5))GO  TO  96 
Fl=1.0/ (1,0-Y (J) ) 

JK=J-1 

F2-1.0/  (1,0-Y (JK) ) 

FB=1.0/(1.0-B) 

FRAC= (F3-F1)/ (F2-F1) 

Jl=L 

DO  99  J2=l,8 
DO  99  J3=l, 5 

99  E(J1,J2,J3)=E(J1,J2,J3)+FRAC*(C(J1,J2,J3)-E(J1,J2,J3) ) 
R(1)=B 

S= ( 1 o  0-B ) /4 . 0 
DO  33  JI=2, 5 
33  R(JI)=R(JI-1)  +  S 
IF (KK-4)  35,100,35 
35  IF (KK-7 )  108,100,108 

96  S=(1.0-Y(J) )/4.0 
R(1)=Y(J) 

DO  97  JI=2, 5 

97  R(JI)=R(JI-1)  +  S 
R(1)=B 

IF (KK-4)  108,100,100 

100  IF (D (L, 1, 1) cLE.0.01*D (L, 1, 5) )GO  TO  102 
J1=L 

DO  101  J2=l, 8 
DO  101  J3=l, 5 

101  F(J1,J2,J3,  )=F  (Jl,  J2 ,  J3)+FRAC*  (D(J1,J2,J3)-F(J1,J2,J3)  ) 

102  IF(E(L,1,1).GE.0.01*E(L,1,5) )  GO  TO  103 

FRAC= (SQRT (A) -SQRT (X (1-1) ) )/ (SQRT (X (I) ) -SQRT (X (1-1) ) ) 

GO  TO  104 
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103  FRAC= (A-X(I-l) )/ (X(I)-X(I-l) ) 

104  Jl=L 

DO  105  J2=l,8 
DO  105  J3=l,  5 

105  E(Jl,J2,J3)=F(Jl,J2,J3)+FRAC*(F(Jl#J2,J3)-E(Jl,J2,J3)) 

108  RTIM=TM*Z(L) 

PRINT  502,Z(L),RTIM 
M=1 

DO109  IL=1,4 
Ml=Mfl 

PRINT  503, W(M)  ,  (E  (L,M,K)  ,K=1,5)  ,W(M1)  ,  (E  (L,M1,K)  ,K=1,5) 
DOlll  K=l,5 

TEMPI  (K)=E  (L,M,K)  *DEMf  TI 
111  TEMP2  (K)=E  (L,M1,K)  *DEMfTI 

PRINT  505, (TEMPI (K) ,K=1, 5) , (TEMP2 (K) ,K=1, 5) 

PRINT  504, <R(K) ,R=1,5) , (R(K) ,K=1,5) 

M=Mf  2 

109  CONTINUE 
11=11+1 

IF(II-Nl)  110,110,115 

110  CONTINUE 
115  REWIND  1 

GO  TO  1 
120  CONTINUE 
STOP 
END 
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PROGRAM  PLT ID (IN PUT -OUTPUT) 

DIMENSION  X(12),Z(25),TIME(25),B(5),TEHP(5),C(25,5),D(25,5), 
1E(25,5),F(25,5),G(25,5),P(25,5),R(25,5),S(25.5),T(25,5) 

2U (25,5)  ,V (25,5)  ,W  (25,  5)  ,Y  (25# 5)  * 

DATA  3i/.01,.02,.05,.10,.20,.50,1.0,2.0,5.,0,10.0,20.0,5p.0/  , 

DATA  Z/0.0,0.5,0.75,1.0,1.2,i.4,1.6^1.8,2.0,2.2,2«5,2.8,3#25,4.0,5 
+.0,6. 0,7. 0,8.0,9*0,10.0, 11. 0,12. 0,14. 0,16. 0,20.0/  1  , 

READ  103, ((E(Kl,K2),Kl=l,25),K2=l,5)  . 

READ  103,  ( (F(K1,K2),K1=1,25),K2=1,5) 

RE^D  103,  (  (G  (Kl,K2)  ,K1=1,'25)  ,K2=1,5) 

READ  103, ( (H(K1,K2),K1=1,25),K2=1,5) 

READ  103,  ( (P(Kl,K2)  ,K1=>1,25)  ,K2=1,5) 

READ  103, ( (R(Kl,K2) ,Kl=l,25) ,K2=1>5)  , 

READ  103, ((S(Kl,K2),Kl*l,25),K2=l,5)  • 

READ  103,  ((T(Kl,K2),Kl=l,25),K2i=l,5) 

READ  103, ((U(K1,K2),K1=1,25),K2=1,5) 

READ  103, ( (V(Kl,K2),Kl=l,25),K2=l,5)  1 

READ  103, ((W(Kl,K2),Kl=l,25),K2=l,5) 

READ  103, ( (Y(Kl,K2),Kl=l,25),K2=l,5)  ! 

1  READ  101 ,  XK,  RHO  ,  CP,  TM,  THICK,  TI  ,N 

101  FORMAT (6E10. 4, 110)  1 

103  FORMAT (10F8.5/10F8.5/5F8. 5)  ’ 

IF(N.EQo0)GO  TO  2 

IF (XK)  100,100,99  1  , 

99  READ  102, QT, ABSR,  (TIME  (M),M=1,N) 

102  FORMAT (2E10.4, 12F5«2/13F5«2) 

GO  TO  5  '  ’ 

2  IF (XK) 100 *  100.7 

7  READ  102, QT, ABSR  ’  i 

DO  4  M=l, 25 

4  TIME  (M)=Z  (M) 

5  ALPHA=XK/ (RH0*CP)  '  , 

A=ALPHA*TM/  (THICK*THICK) 

PRINT  401, A, THICK 

401  FORMAT (1H1,36H  BASIC  DIMENSIONLESS  PARAMETER  BPLT= , F8 • 5 , 32X, 16HPLA 
1TE  THICKNESS* , F6 • 3 , llHCENTIMETERS ) 

Pn  INT  402  XK  RHO  CP 

402  FORMAT ( 1H* , 28H  MATERIAL  THERMAL  PROPERTIES , 4X, 2 lHTHERMAL  CONDUCTIV 
1 ITY= , F6 . 3 , 12HCAL/SEC-CM-C , 4X, 8HDEN!3ITY= , F5 . 2 , 7HGMS/CM3 , 4X, 14HSPECI 
2FIC  HEAT= , F4 . 3 , 8HCAL/GM-C ) 

PRINT  403,TM,QT 

403  FORMAT (1H,31H  WEAPON  THERMAL  CHARACTERISTICS, 12?, 2 3HTIME  TO  ,MAX  I 

1RRADIANCE= ,  F7 . 4 , 7HSECONDS ,  9X,  2  3IITOTAL  RADIANT  EXPOSURE* ,  F8 . 2 , 7HCAL 
2/CM2 )  i 

HM=QT/(2.6*TM) 

DEM=ABSR*HM*THICK/XK  ' 

B(l)=0o0 

DO  6  M=2 , 5  i 

M1=M-1  i 

6  B(M)=B(M)+THICK/4,0 

Figure  3.  Computer  Program  to  Obtain  the  Temperature  Distribution, 
History  in  Plates 
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LT«1 

DO  9  J=l,  12 
I=J 

IF(  A-X  ( J)  )  10,10,9 
.  9  CONTINUE 

1*13 

10  GO  TO(51,$2, 53,54,55,56,57, 58,59,70,71,72, 73)1 

51  DO  61  16=1,25 
DO  61  N=l,  5 

61  C  (H,N)=E  (M,N) 

LL=0 

GO  TO  90 

52  DO  62  16=1,25 
DO  62  N=l, 5 
C(M,N)=E(M,N) 

62  D(M,N)=F(M,N) 

GO  TO  90 

53  DO  63  K=l,25 
«  DO  63  N==l,  5! 

C(M,N)*F(M,N) 

63  D(M,N)=G(M,N) 

GO  TO  90 

54  DO  64  16=1,25 
DO  64  .N=l,5 
C(M,N)=G(M,N) 

64  D<M,N)=H(M,N) 

GO  TO  90 

55  DO  65  16=1,25 
DO  65  11=1, 5 

C  (M,N)=H  (M,N) 

65  D(M,N)=P(M,N) 

GO  TO  90 

56  DO  66  M-l, 25 
DO  66  N=l,5 
C(Mi,N)=P(M,N) 

66  D(M,N)=R(M,N) 

GO  TO  90 

57  DO  67  16=1,25 
DO  67  N=l, 5 

,  C (M,N)=R (M,N) 

67  D(M,N)=S(M,N) 

GO  TO  90 

58  DO  68  N=l, 25 

DO  68  N=l,  5 
C(M,N)=S(M,N)  . 

68  D(M,N)=T(M,N) 

GO  TO  90 

59  DO  69  M=l, 25 
DO  69  N=l, 5 

C (M,N)=T (M,N) 

69  E}(M,N)=U(M,N) 

GO  TO  90 

7p  DO  80  M=l;25 
DO  80  N=l, 5 
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C(M,H)=U(M,H) 

80  D(M,H)=V(M,N) 

GO  TO  90 

71  DO  81  M=l,25 
DO  81  *=  1,5 
C(M,N)=V(M,H) 

81  D(M,H)=1«(IM0 
GO  TO  90 

72  DO  82  11=1,25 
DO  82  »=1,5 
C(M,N)=W(M,H) 

82  D  (M,H)=Y  (M,H) 

GO  TO  90 

73  DO  83  11=1,25 
DO  83  11=1.5 

83  C  CM,HJ=Y (M,N) 

LL=G 

90  11=1 

DO  50  1=1,25 

IF  (TIME  (II)— Z  (L)  )  25,25,50 
25  IF(LL)  35,35,30 
30  IF (C (L, 1) ,GE.o.ol*C (L, 5) )GO  TO  32 

FRAC=(SQRT(A)-SQRT(X(I-1)  )  )/(SQRT(X(I)  )-SQRT(X(I-l)  )  ) 

GO  TO  33 

32  FRAC=(A-X{I-1))/(X(I)-X(I-1)) 

33  DO  34  J=l,5 

34  C  (L, J)=C  (L,  JJ+FRAC*  (D  (L, J)  — C  (L, J) 

35  RTIM=TM*Z  (L) 

PRINT  202, Z (L) ,RTIM 

202  FORMAT (1H0,2X,42H  TEMPERATURES  FOR  A  DIMENSIONLESS  TIME  OF  ,F5.2, 
130X,6H  TXME=,F6.3,7HSECQNDS) 

DO  40  J=l,  5 

40  TEMP  (J)=C  (L, J)  ♦DEM'^TI 
PRINT  203, (B(J),J=1,5) 

203  FORMAT (1H0,8X,12H  DEPTH(CM.)=,5F9.2) 

PRINT  204, (TEMP (J) , J=l, 5) 

204  FORMAT (1H  ,3X,17H  TEMPERATURE (  C)=,5F9.2) 

50  CONTINUE 

GO  TO  1 
100  CONTINUE 
STOP 
END 
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program  Karen.  (atvr.aumr) 

DIMENSION  X(9)  ,Y(5) ,W(12) ,Z(24) ,R{5)  ,C (24. 12. 5)D(24. 12. 5) .1 (24. 12 
1.  5)  .F  (24. 12.  5)  .TIME  (24)  .TEMPI  (5)  .TEMP2  (5) 

DATE  X/.010. .025. .050. .100/ 

DATA  Y/0. 95, 0.90, 0.80/ 

DATA  W/0.. 30.. 60.. 90.. 120.. 150.. 180.. 210.. 240. .270.. 300.. 330./ 

DATA  Z/0. 5. 0.75. 1.0. 1.2. 1.4. 1.6. I.8.2.O. 2. 2. 2.5. 2.8. 3.25. 4.00. 5.0. 
-1-6. 0.7. 0.8.0.  9.0. 10.0. 11. 0.12.0. 14. 0.16.0.  20.0/ 

1000  READ  1001.XK.RB0.CP.TM.RI.RO.TI.KK.I.J.H 

1001  FORMAT (7F10.4, 2 (13, 12) ) 

IF(N.EQ.0)O0  TO  2 

11  I(XK)  120.120.1100 
1100  READ  1002. OT. ABSR.  (TIKE (M).Mb1.B) 

1002  FORMAT  (2E10.4.12F5.2/13F5.2) 

502  FORMAT (1H0.2X.56B  DIMENSIONLESS  TEMFERATORES  FOR  A  DIMENSIONLESS  T 
1IME  OF  .F5.2.30X.6H  TXMB«.F6.3,7BSS0GHDS) 

503  FORMAT (1H0.11H  CYL  ABGLfc*,F4.0,  3HDEG,3X,5HT*  =,5F7.4.5X.  10HCYL  AS 
1GLE=,F4.0. 3HDEG.2X.5HT*  -.5F7.4) 

504  FORMAT  (1H  ,  2lX,  5HR/RO=,  5F7  .4.  24X,  5HR/RO=,  5F7  .4) 

505  FORMAT (1H  ,2lX, 5HT (C)=. 5F7.1.24X, 5HT(C)=.5F7.1) 

GO  TO  5 

2  IF(XK)  120,120,7 
7  READ  1002 , QT, ABSR 
DO  4  M=l,24 

4  TIME  (M)  =A  (M) 

Nl=24 

5  ALPHA=XK/  (RHO*CP) 

A=ALPHA*TM/ (RO*RO) 

B=RI/RO 

PRIST  401, A, B 

401  FORMAT91H1.31H,  BASIC  DIMENSIONLESS  PARAMETERS, 20X,5HBCYL=,F7. 6, 
15X,  6HRI/R0=,  F4 . 3) 

PRINT  402.XK.RHO.CP 

402  FORMAT  (1H  ,28H  MATERIAL  THERMAL  PROPERTIES, 4X,  2 1HTHERMAL  CONDUCTIV 
1ITY=,F6.3, 1 2HCAI/  SEC-CM-C  ,  4X,  8HDENSITY=  ,  F5 . 2 , 7HGMS/CM3 , 4X,  14HSPECI 
2FIC  HEAT= , F4 . 3 , 8HCAL/GM-C ) 

PRINT  403.TM.QT 

403  FORMAT  (1H.31H  WEAPON  THERMAL  CHARACTERISTICS,  10X.23HTIME  TO  MAX  I 
1RRADIANCE=  , F7 . 4 , 7HSEC0FDS , 9X, 2 3HT0TAL  RADIANT  EXPOSURE^, F8.2, 7HCAL 
2/CM2) 

PRINT  404,RI,RO,TI 

404  FORMAT (1H, 27X, 13HINNER  RADIUS= , F6 . 3 , llHCENTIMETERS , 4X, 1 3H0UTER  RA 
1DIUS= , F6 . 3 , 1 lHCENT IMETERS , 4X, 2  OHINITIAL  TEMPERATURE^ , F5 . 1 , 5HDEG-C ) 

IZ=0 

HM=QT/ (2 . 6*TM) 

DEM-ABSR*HM*RO/XK 

R(1)=B 

S=(1.0-B)/4o0 
DO  3  K=2,5 
Kl=K-l 

Figure  4„  Computer  Program  to  Obtain  the  Temperature  Distribution 
History  in  Circular  Cylinders  Rotating  at  6  rpm. 
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3  R{K)=R(K)+S 

45  GO  TO  (50, 60,70,60, 60, 70,€O)XX 

50  READ  1.  (((EiKl.X2.K3)  ,10*1,24)  ,K2=1,12)  ,K3=1.5) 

1  FORMAT (8F10. 7) 

GO  TO  83 
60  DO  21  K3»l,5 
DO  21  K2=l,12 

21  READ  1,  (C(Kl,K2,K3),Kl=l,24) 

DO  22  K3»l,5 
DO  22  K2=1,12 

32  READ  1,  (D(Kl,K2,K3),Kl=l,24) 

IF (KK— 4)  88,80,80 

70  READ  1,  (  (  (C  (Kl,K2,K3)  ,Kl=l, 24)  ,K2=1, 12)  ,K3=1, 5) 

READ  1,  (((E(K1,K2,K3),K1^1,24),K2*1,12).K3*1.5) 

GO  TO  88 

80  READ  l,(((E(Kl,K2,K3),Kl*il,24),K2*l.l2),K3*l,5) 

READ  1,  (  (  (F  (K1,K2,K3)  ,Kl=l, 24)  ,K2=1, 12)  ,K3=1, 5) 

88  11=1 

DO  110  L=l,24 

IF  (TIME  (XI)— Z  (L)  )  90,90,110 

90  GO  TO  {108,  91,  98,  98,  98.  98) KK 

91  IF (C (L, 1, 1) ,GE.0,01*C (L, 1, 5) )G0  TO  92 

FRAC= (SQRT (A)-SQRT (X(I-l) )  )/  (SQRT (X(I) )-SQRT (X(I-l) ) ) 
GC  TO  93 

92  FRAC=(A-X(I-1))/(X(I)-X(I-1)) 

93  J1=L 

DO  94  J2=l, 12 
DO  94  J3=l,5 

94  E(J1,J2,J3)=C(J1.J2,J3)+FRAC*(D(J1,J2,J3)-C(J1,J2,J3)) 
GO  TO  108 

98  IF(C  (L,  1,1)  „LE.0.01*C(L,  1, 5)  )G0  TO  96 
Fl=l,0/  (1,0— Y  (J)  ) 

JK=J-1 

F2=1.0/(1.0-Y(JK)) 

FB=lo0/ (1,0— B) 

FRAC= (FB-Fl)/ (F2-F1) 

J1=L 

DO  99  J2=l, 12 
DO  99  J3=l, 5 

99  E  (Jl,  J2,  J3)=E (Jl, J2, J3)+FRAC* {C (Jl, J2, J3)-E (Jl,  J2,  J3) 
R(1)=B 

S= (l„0-B)/4o0 
DO  33  Jl=2 , 5 

33  R  (JI)=R  (JI-1)+S 

IF (KK-4)  35,100,35 

35  IF (KK-5)  36,100,36 

36  IF (KK-7 )  108,100,108 

96  S=  ( 1 . 0-Y ( J ) ) /4 , 0 
R(1)=Y{J) 

DO  97  JI=2, 5 

97  R(JI)=R(JI-1)+S 
R(1)=B 

IF (KK-4)  108,100,100 
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100  XFCKK-61  41,108.41 

41  XF (0(1^  1,1)  •LB«0I,C1*D(I<,1,5))G0  TO  102 
Jl mL 

DO  101  J2»lf12 
DO  101  J3»l,5 

101  F(J1,  J2,  J3)»F(Jl,  J2,J3)+FmC*  (D(Jl,  J2,  J3)-F  (Jl,  J2,  J3)  ) 

102  IF(E  (L,  1,1)  ,GE,0,01*E  (£>,1,5)  GO  TO  103 
FRAC=(SQI»(A)-S<X*TC*(I-1)  )  )/  (SORT  (X  (I)  )-SQRT  (X(X-l)  )  ) 
GO  TO  104 

103  FRAC*(A-X(I~1)  )/(X(I)-X(I-l)) 

104  Jl*L 

DO  105  J2=i,l2 
DO  105  J3=l,5 

105  F  (Jl, J2, J3)=E  (Jl, J2, J3J+FRAC*  (F  (Jl, J2, J3J-E  (Jl, J2, J3)  ) 

108  RTI**TH*Z(L) 

FRIST  502,Z(L),STXM 

Hel 

D0109  UE=1,6 
lZ=IZfl 

PRIST  503,W(H),E(L,M.K),K=1,5),W{M1),  (E(L,K1,K)  ,K=1.5) 
DOlll  K=l,5 

TEMPI  (K)=E  (L,M,K)  *DEHfTI 
111  TEMP2  (K)=E  (L,Ml,K)  *DEK+TI 

PRIST  505,  (TEMPI (K)  ,K»1,5)  ,  (TEMP2  (K)  ,K=1,5) 

PRIST  504,  (R(K)  ,R=1,5)  ,  (R(K),K=1,5) 

M=M+2 

109  CONTINUE 
IF(IZ-12)  131,130,130 

130  IZ=0 
PRIST  509 

509  FORMAT (lHl.//) 

131  CONTINUE 
11=11+1 

IF(II-Nl)  110,110,115 

110  CONTINUE 
115  GO  TO  1000 
120  CONTINUE 

STOP 

END 
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PROGRAM  TOWEL  (MPOT.OOTPOT) 

DDCKSKS  T(30)  .Q(30)  ,TEM?(30)  .TIME  (30) 

!MTA  l/0#0*0#5|0«75|l«0fl»2»l»4»l»6#l#8#2»0»2»2#2»5#2»8*3»25»4#0, 
+5. 0,8, 0,7. 0,8.0, 9.0, 10. 0,12. 5. 15. 0,20.0,  30. 0,40. 0,50.0/ 

DATA  0/0.0.  .037.  .109.  .203.  .279.  .347.  .400.  .422.  .477.  .507.  .545.  .577 
♦  .017.  .468,  .714.  .748.  .773.  .793.  .810.  .824.  .849.  .874.  .902.  .934.  .952.  . 
*964/ 

1  REID  100.  SBO.CP.TIi.  THICK.  TI.QT.AE68 
100  FORMAT  (7P10. 4) 

IF(KHO)  75.75.2 

2  DE8fc=RH0*CP*THICK 
DO  20  R=l,26 

20  T«(X)*T(K)*TM 
PRINT  200 

200  FORMAT  (1H),  74H  THERMALLY-THIN  ELEMENTS  AVERAGE  TEMPERATURE  HIS 
1TORY  BY  ENERGY  BALANCE) 

PRINT  201, RBO, CP, THICK 

201  FORMAT (1H0.11H  DENSITY- , F6 . 3 , 7HGMS/CM3 , 14H  SPEC.  HEAT=,F5.3, 
112HCAL/CM-SEC-D,  15H  MALL  THICK.=,F6.3,  3HCM.) 

PRINT  403,TM, QT 

403  FORMAT (1h',31H  WEAPON  THERMAL  CHARACTERISTICS,  3X.23HTIME  TO  MAX  I 
lRRADIANCE-  ,  F6  .  4  , 7HSECCBDS ,  4X,  2 3HTOTAL  EXPOSURE-, F6. 2, 7BCAL 
2/CM2) 

DO  10  K=l,26 

10  TEMP  (K)=Q  (K)  *QT/DEMfTI 
PRINT  204.  (T  (L)  ,I#=1, 14) 

204  FORMAT (1H0.12H  T/TMAX=, 14F7.2) 

PRINT  205,  (TIME (L),  1^1, 14) 

205  FORMAT (1H,  12H  TIME (SEC)=, 14F7.3) 

PRINT  206, (TEMP(L) ,L=1,14) 

206  FORMAT (1H  ,12H  TEMP(  C)=,14F7.1) 

PRINT  207, (T(L),L=15,26) 

207  FORMAT (1E0.16X,10H  T/TKAX-, 12F7 .2) 

PRINT  208, (TIME (L) ,L=15,26) 

208  FORMAT (1H  ,16X,10HTIME(SEC)=,12F7.3) 

PRINT  209,  (TEMP(L),L=15,26) 

209  FORMAT (1H,  16X.10H  TEMP (  C)=,12F7.1) 

GO  TO  1 

75  CONTINUE 
STOP 
END 


Figure  5C  Computer  Program  to  Obtain  the  Temperature  History  in 
Thermally -Thin  Elements 
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BASIC  DIMENSIONLESS  PARAMETER,  BELT 
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0  2.5  5.0  7.5  10.0  12.5  15.0 

DIMENSIONLESS  TIME,  »* 


FIG.  6  THERMAlLY-THICK  AND  THERMALLY-THIN  REGIONS  OF  A  PLATE 
EXPOSED  TO  THE  THERMAL  RADIATION  PULSE  OF  A  NUCLEAR 
WEAPON. 
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FIG.  7  THERMAIN-THICK  AND  THERMAIN-THIN  REGIONS  OF  A  CIRCULAR  CYLINDER 
EXPOSED  TO  THE  THERMAL  RADIATION  PULSE  OF  A  NUCLEAR  WEAPON 
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FIG.  7 (CONTINUED) 
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NUMERICAL  CALCULATION 
INTERPOLATION  RULE  (EQUATION  (17)  ) 
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FIG.  8  COMPARISON  OF  THE  SURFACE  TEMPERATURE  ON  A  STATIONARY  CIRCULAR  CYLINDER 
AS  GENERATED  BY  A  NUMERICAL  CALCULATION  AND  BY  THE  INTERPOLATION  RULES. 
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APPENDIX  A 

HEAT  TRANSFER  EQUATIONS  AND  BOUNDARY  CONDITIONS 


The  heat  conduction  equation  is  the  governing  partial 
differential  equation  which  determines  the  temperature  distribution 
history  in  the  structural  elements  of  this  report.  These  elements 
are  heated  by  the  nuclear  weapon  thermal  radiation  pulse  only  and  no 
heat  is  removed  at  any  surface.  Also,  these  structural  elements  are 
assumed  to  be  constructed  of  isotropic,  homogeneous  materials  whose 
thermophysical  properties  are  constant.  Following  is  a  mathematical 
statement  of  the  thermal  radiation  heating  problem  for  each  of  the 
four  structural  elements,  consisting  of  the  heat  conduction  equation 
and  appropriate  boundary  conditions. 


(1)  Plates 
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(2)  Circular  Cylinders 
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(3)  Rotating  Circular  Cylinders 
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Bote  that  afeti/S  radians/sec  (6  rpm)  only  in  this  report. 
(4)  Thermal ly-Th in  Elements 
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Figure  (A-l)  is  a  schematic  diagram  of  each  of  the  four 
structural  elements  of  this  report.  Also  shown  are  the  positions 
on  the  element  for  which  temperature  has  been  calculated  by 
solving  the  appropriate  heat  conduction  equation  described  above. 
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APPENDIX  B 

DEFINITION  OF  VARIABLES  REQUIRED  TO  RUN  THE  COMPUTER  PROGRAMS 


*  i 

The  variables  used  in  the  four  computer  programs  of  this  report 
have  common  definitions  and  use. common  systems. of  units.  Two  - 
possible  systems  of  units  may  be  used.  If  temperatures  in  degrees 
Centigrade  are  desired,  all  variables  are  given,  in  c.g.s.  units  and 
calories.  If  degrees  Fahrenheit  are  wanted,  the  corresponding  input 
variable  units  are  inches,  seconds  and  BTU.  The  following  is  a  list  1 
of  definitions  for  all  of  the  variables  required'  to  run  any  of  the 
four  computer  programs. 

I 

Floating  Point  Variables 

•  1 

Absorptance  of  target  material  ; 

Material  specific  heat ; 

Material  thermal  conductivity 
Inner  radius  of  the  circular  cylinders 
Material  density  , 

Outer  radius  of  the  circular  cylinders 
Initial  temperature 

Dimensionless  tiime  values  for  which  the  temperature  history  ’ 
is  to  be  calculated,  (chosen  from  the  allowable  values  of 
Table  1)  1  ■ 

Plate  thickness 

Time  at  which  maximum  thermal  irradiance  occurs  1 
Total  radiant  exposure  per  unit  area 

.  i 

Fixed  Point  Variables  1  1 

-  ,  i  1 

I  Number  indicating  the  value  of  BCYL  for,  rotating  cylinders 

(see  description  of i program  ROTCYL) 

J  Number  indicating  the  value  of  Ri*  for  rotating  cylinders 

(see  description  of  program  ROTCYL) 1 

KK  Number  indicating  the  manner  of  interpolation  or  extrapolation 

of  circular  cylinder  data  (see  description  of  program  ROTCYL) 


ABSR 

CP 

XK 

RI 

RHO 

RO 

TI 

TIME 


THICK 

TM 

QT 


I 


*  The  user  is  cautioned  that  the  headings  on  the  computer  outputs 
are  always  in  c.g.s.  units,  calories  and  degrees  centigrade. 
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APPENDIX  C 

EXAMPLES  FOR  THE  SET-UP  AND  USE  OF  THE  COMPUTER  PROGRAMS 


“jOiis  section  presents  examples  showing  how  each  of  the  computer 
programs  of  this  report  are  set-up  and  run.  The  output  of  each 
■program  will  also  be  described  and  part  of  the  solutions  presented 
in  graphical  form.  A  single  target  element  will  be  analyzed  by  each 
program  in  turn.  It  was  chosen  such  that  the  results  of  each  program 
can  be  compared  to  show  when  these  solutions  are  in  agreement. 

1  Figures  6  and  7  which  indicate  where  heat  conduction  effects  become 
important  will  be  referred  to,  to  explain  why  these  programs  can  yield 
the  same  result.  The  common  target  element  is  composed  of  aluminum 
which  has  the  following  thermal  properties:  Thermal  conductivity  = 
0.374  cal/sec-cm-°C,  Density  =  2.71  gms/cm3,  and  Specific  Heat  = 

0.23  cal/gm-°C.  (Thermal  Diffusivity  =  0.60  cm2/sec) a  This  element 
is  normally  exposed  to  a  nuclear  weapon  detonation  such  that  the 
time  to  maximum  irradiance  is  1.0  second  and  is  at  a  separation 
distance  from  the  explosion  such  that  the  total  radiant  exposure  is 
100  cal/cm2.  The  element  initial  temperature  is  20°C  and  its 
absorptance  is  1.0. 

(1)  Temperature  Distribution  History  in  a  Stationary  Circular 
Cylinder  (Program  CYLDAT) 

A  cylindrical  target  element  of  outer  radius  6.33  centimeters 
and  inner  radius  5.38  centimeters  is  exposed  to  the  nuclear  weapon 
detonation  described  above.  Two  data  input  cards  plus  the  magnetic 
tape  containing  the  dimensionless  temperature  data  are  required  to 
run  this  program.  The  following  is  a  table  of  the  numbers  required 
and  their  location  on  the  input  cards.  The  corresponding  variables 
were  specified  in  the  Description  and  Use  of  Computer  Programs  section 
and  were  defined  in  Appendix  A. 

,  Card  Column  Number 

10  20  30  40  50  60  70  80 

Card  1  0.374  2.71  0.230  1.0  5.38  6.33  20.0  0 

Card  2  100.0  1.0 

i 

For  this  input,  the  program  computes  and  prints  out  the  values  of 
BCYL  <*nd  Ri*  (0.015  and  0.85,  respectively).  The  values  of  the 
material  thermophysical  properties  (k,  p,  and  cp),  the  weapon  thermal 
pulse  characteristics  (tm  and  QT)  and  the  cylinder  radii  are  also 
printed.  Subsequently,  the  temperature  history  is  retrieved  and 
printed  in  dimensional  and  dimensionless  form  for  the  24  selected 
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values  of  tine,  the  eight  selected  values  of  cylinder  angle  and  the 
five  allowable  values  of  cylinder  radius  (see  Table  1).  Figure  B-l 
is  a  plot  of  part  of  the  results  generated  by  the  program  in 
solving  this  problem. 

(2)  Temperature  History  in  a  Plate  (Program  PI/TlD ) 

A  target  element  consisting  of  a  plate  0.95  centimeters  thick 
is  exposed  to  the  sane  nuclear  weapon  detestation  as  in  the  previous 
example.  Note  that  this  plate  has  the  same  thickness  as  the  stationary 
circular  cylinder.  Two  data  input  cards  are  required  to  run  this 
program.  The  following  is  a  table  of  the  numbers  and  their  location 
on  the  data  input  cards.  The  corresponding  variables  were  specified 
in  the  Description  and  Use  of  Computer  Programs  section  and  were 
defined  in  Appendix  A. 


Card  Column  Number 

10  20  30  40  50  60  70 

Card  1  0.374  2.71  0.230  1.0  0.95  20.0  0 

Card  2  100.0  1.0 

For  this  input,  the  program  prints  the  value  of  BPLT  and  the  plate 
thickness.  The  values  of  the  thermophysical  properties  (k,  p  and  cp  ) 
and  the  thermal  pulse  characteristics  (tm  and  QT)  are  also  printed. 
Subsequently  the  temperature  history  is  retrieved  and  printed  for 
the  25  selected  values  of  time  and  the  five  selected  values 
of  distance  into  the  plate.  Figure  B-2  is  a  plot  of  part  of  the 
results  generated  by  the  program  in  solving  this  problem. 

A  comparison  of  figures  B-l  and  B-2  shows  that  the  temperature 
history  of  the  circular  cylinder  at  its  most  forward  point  ($=0°) 
is  similar  to  that  of  the  plate  with  identical  thickness.  This 
result  could  have  been  predicted  by  referring  to  Figure  7.  That  is, 
for  the  values  of  BCYL  and  Ri*  of  this  problem  (.015  and  0.85), 

Figure  7  shows  that  angular  conduction  is  negligible  for  t*  values 
up  to  2.0.  A  t*  of  2.0  corresponds  to  2.0  seconds  of  real  time  for 
this  problem.  When  angular  (circumferential)  heat  conduction  is  small, 
the  most  forward  cylinder  point  and  the  plate  will  have  similar 
temperature  histories.  As  angular  conduction  increases  v’ith  time, 
more  heat  is  conducted  away  from  the  most  forward  point  and  its 
temperature  history  departs  from  that  of  a  plate. 

(3)  Temperature  Distribution  History  in  a  Rotating  Circular  Cylinder 
(Program  ROTCYL) 

The  cylindrical  target  element  of  example  (1)  is  rotating  at  six 
revolutions  per  minute  and  is  exposed  to  the  same  nuclear  weapon 
detonation.  Before  this  program  can  be  run  it  is  first  necessary 
to  compute  the  values  of  the  parameters  BCYL  and  Ri*.  Values  of 
0.015  and  0.85  are  found  for  these  variables.  Next,  these  values  must 
be  inspected  to  determine  if  interpolation  or  extrapolation  of  the 
dimensionless  temperature  data  is  required.  An  inspection  shows 
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that  interpolation  cm  both  BCYL  and  Ri*  is  required  (see  Table  1 
for  selected  values).  Furthermore,  KK,the  number  indicating  the 
manner  of  the  interpolation  is  four  (see  Table  under  Description 
and  Use  of  the  Computer  Programs).  Finally,  the  dimensionless 
temperature  decks  must  be  selected  and  the  numbers  I  and  J  indicating 
these  decks  specified.  Z  indicates  the  nearest  larger  value 
of  BCYL  and  j  the  nearest  smaller  value  of  Ri*  hence  I  is  two  and  J  is 
three  for  this  case  (allowable  values  of  I  and  J  are  specified  in 
Description  and  Use  of  the  Computer  Programs).  Some  double 
interpolation  is  needed  in  this  problem,  four  dimensionless 
temperature  data  decks  are  input  in  addition  to  the  two  data  input 
cards.  The  following  is  a  table  of  the  numbers  required  and  their 
location  on  the  input  cards.  The  corresponding  variables  were 
specified  in  the  Description  and  Use  of  Computer  Programs  section 
and  were  defined  in  Appendix  A. 

Card  Column  Humber 

10  20  30  40  50  60  70  73  75  78  80 

Card  1  0.374  2.71  0.230  1.0  5.38  6.33  20.0  4230 

Card  2  100.0  1.0 

Dimensionless  Temperature  Data  Deck  (1=1,  J=3) 

Dxmensionless  Temperature  Data  Deck  (1=2,  J=3) 

Dimensionless  Temperature  Data  Deck  (1=1,  J=2) 

Dimensionless  Temperature  Data  Deck  (1=2,  J=2) 

For  this  input,  the  program  prints  the  value  of  BCYL  and  Ri*.  The 
values  of  the  material  thermophysical  properties  (k,  p ,  and  cp),  the 
weapon  thermal  characteristics  (tm  and  QT)  and  die  cylinder  radii 
are  also  printed.  Subsequently,  the  temperature  distribution  history 
is  retrieved  and  printed  in  dimensional  and  dimensionless  form  for 
the  24  selected  values  of  time,  the  12  allowable  values  of 
cylinder  angle  and  the  five  selected  values  of  cylinder  radius 
(see  Table  1) .  Figure  B-3  is  a  plot  of  part  of  the  results  generated 
by  the  program  in  solving  this  problem. 

Figures  B-l  and  B-3  can  be  compared  to  show  the  effect  of 
cylinder  rotation  on  temperature  distribution  history.  It  is 
seen  that  rotation  causes  great  changes  by  (1)  lowering  the  maximum 
temperatures  (2)  decreasing  the  temperature  gradients  in  the  circumfer¬ 
ential  direction  and  (3)  causing  the  maximum  temperatures  to  occur  at 
a  cylinder  angle  of  30°  instead  of  at  the  most  forward  point.  A 
comparison  of  these  figures  with  Figure  (B-2)  shows  that  the 
temperature  history  of  the  most  forward  point  of  cylinders  and  the 
plate  are  similar  for  only  the  first  second  of  heating.  Thus,  the 
region  of  negligible  angular  conduction  shown  in  Figure  6  for  a 
stationary  cylinder  has  been  reduced  to  approximately  one  second  by 
the  6  rpm  rotation.  Also,  it  can  be  deduced  from  the  figures  that 
rotation  causes  a  much  greater  departure  from  the  temperature  history  for 
no  angular  conduction  than  a  movement  into  the  angular  conduction 
region. 
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(4)  Temperature  History  in  a  Thenaally-Thin  E lament  (Program  TBHEL) 

A  target  element  0.95  centimeters  thick  is  exposed  to  the 
nuclear  weapon  detonation  of  the  previous  examples.  Bote  that  this 
e lament  can  represent  either  the  plate  or  cylindrical  element  of 
examples  (1)  and  (2).  Only  one  data  input  card  is  required  to  run 
this  program.  The  following  is  a  table  of  numbers  and  their 
location  on  the  data  input  card.  The  corresponding  variables  were 
specified  in  the  Description  and  Use  of  Computer  Programs  section 
and  were  defined  in  Appendix  a. 

Card  Column  Humber 

10  20  30  40  50  60  70 

Card  1  2.71  0.230  1.0  0.95  20.0  100.0  1.0 

For  this  input,  the  program  prints  the  values  of  the  material 
thermophysical  properties  ( 0,  and  c_) ,  the  weapon  thermal  pulse 
characteristics  (tm  and  QT)  and  the*element  thickness.  Subsequently , 
the  average  temperature  history  is  confuted  and  printed  for  the 
26  selected  values  of  time  (see  Table  1) .  Figure  B-4  is  a  plot  of 
the  results  generated  by  this  program. 

A  comparison  of  Figures  B— 2  and  B-4  shows  that  the  temperature 
history  of  the  plate  and  the  thin  element  are  very  similar  for 
large  values  of  time.  This  result  could  have  been  predicted  by 
Figure  6.  That  is.  Figure  6  predicts  that  the  plate  is  thermally- thin 
for  dimensionless  time  values  greater  than  3„5  (greater  than  3.5 
seconds  in  these  examples) .  A  comparison  of  Figures  B-i  and  B-4 
shows  that  the  temperature  history  of  the  most  forward  point  on  the 
circular  cylinder  and  the  thin  element  are  similar  for  some 
intermediate  values  of  time.  This  region  of  agreement  could  have 
been  deduced  from  Figure  7  as  the  region  where  the  cylinder  is 
both  thermally- thin  and  where  angular  conduction  is  negligible. 
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FIG.  C-l  EXAMPLE  OF  THE  TEMPERATURE  DISTRIBUTION  HISTORY  IN  A  STATIONARY 
CIRCULAR  CYLINDER. 


FIG.  C-2  EXAMPLE  OF  THE  TEMPERATURE  DISTRIBUTION  HISTORY  IN  A  PLATE 


FIG.  C-3  EXAMPLE  OF  THE  TEMPERATURE  DISTRIBUTION  HISTORY  IN  A  CIRCULAR 
CYLINDER  ROTATING  AT  6  RPM. 
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